Ultrasound-induced acoustophoretic motion of microparticles in three dimensions 
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We derive analytical expressions for the three-dimensional (3D) acoustophoretic motion of spher- 
ical microparticles in rectangular microchannels. The motion is generated by the acoustic radi- 
ation force and the acoustic streaming-induced drag force. In contrast to the classical theory of 
Rayleigh streaming in shallow, infinite, parallel-plate channels, our theory does include the effect 
of the microchannel side walls. The resulting predictions agree well with numerics and experimen- 
tal measurements of the acoustophoretic motion of polystyrene spheres with nominal diameters of 
0.537 jim and 5.33 jim. The 3D particle motion was recorded using astigmatism particle tracking 
velocimetry under controlled thermal and acoustic conditions in a long, straight, rectangular mi- 
crochannel actuated in one of its transverse standing ultrasound-wave resonance modes with one 
or two half-wavelengths. The acoustic energy density is calibrated in situ based on measurements 
of the radiation dominated motion of large 5-um-diam particles, allowing for quantitative compar- 
ison between theoretical predictions and measurements of the streaming induced motion of small 
0.5-pm-diam particles. 

PACS numbers: 43.25.Nm, 43.25.Qp, 43.20.Ks, 47.15.-x 



I. INTRODUCTION 

Acoustofiuidics is gaining increasing interest in lab-on- 
a-chip and microfiuidics applications. Techniques based 
on acoustofluidic forces permit to perform a large vari- 
ety of different tasks such as trapping, separation and 
sorting of cells, particle manipulation, and generation of 
fluid motion in a non-intrusive way [1, 2]. Acoustic forces 
allow for non-destructive and label-free particle handling 
based on size, density, and compressibility. Experimen- 
tally, the acoustophoretic motion of particles is driven by 
an ultrasonic standing wave that generates acoustic ra- 
diation forces on the particles and acoustic streaming in 
the fluid, which then exerts a Stokes drag force on the 
particles. Theoretically, such phenomena are described 
by complex, non-linear governing equations sensitive to 
the boundary conditions and are thereby difficult to pre- 
dict. Therefore, the development of analytical and nu- 
merical methods that are able to accurately predict the 
acoustophoretic motion of different particle or cell types 
is currently a major challenge in the design of acoustoflu- 
idic systems [3]. 

To guide and control these theoretical developments, 
precise experimental measurements of the acoustopho- 
retic motion of microparticles are necessary, and particle- 
based velocimetry techniques are among the best meth- 
ods available. The work of Hagsater et al. [4] was one of 
the first to use micro particle image velocimetry (uPIV) 
in resonant microfluidic chips. In their case the measure- 
ments were employed to visualize the resonance modes in 
the microchip, using the radiation-dominated horizontal 
motion of 5-um-diam particles and the associated hori- 



zontal acoustic streaming pattern using 1-um-diam par- 
ticles. Using a similar jiPIV technique, Manncberg et al. 
[5] characterized multiple localized ultrasonic manipula- 
tion functions in a single microchip. Barnkob et al. [6] 
and Koklu et al. [7] also studied acoustophoretic motion 
of large particles (5- and 4-um-diam, respectively), but 
instead used particle tracking velocimetry (PTV) to ob- 
tain particle paths, which were compared with theoret- 
ical results. Later, Augustsson et al. [8] employed both 
PTV and pPIV to make high-accuracy measurements of 
the acoustic energy density as well as the temperature 
and frequency dependence of acoustic resonances in mi- 
crochannels filled with 5-um-diam particles dominated by 
the radiation force. Such approaches have successfully 
been applied to the two-dimensional (2D) motion of par- 
ticles in the optical focal plane in simple geometries and 
resonances. Recently, Dron et al. [9] used defocusing of 
particle images to measure the magnitude of radiation- 
dominated acoustophoretic particle velocities parallel to 
the optical axis in similar simple half- wave resonance sys- 
tems. However, in more complex configurations, or in the 
case of small particles dragged along by acoustic stream- 
ing rolls, more advanced techniques are necessary, that 
are able to resolve three-dimensional (3D) particle posi- 
tions and three-component (3C) motion. Among these 
techniques, those based on jiPIV have issues regarding 
the depth of correlation between adjacent planes [10, 11], 
while classical 3D particle tracking techniques require ei- 
ther stereo-microscopes with tedious calibration proto- 
cols, or fast confocal microscopes with a great loss in 
light intensity due to the use of pinholes [12]. 

In this work, an analytical and experimental anal- 
ysis is presented with the aim to improve the under- 
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standing of the full 3D character of ultrasound-induced 
acoustophoretic motion of microparticles. In particu- 
lar, analytical expressions for this motion are obtained 
by extending the classical results for Rayleigh streaming 
in shallow parallel-plate channels to also cover rectan- 
gular channels of experimental relevance. The analyti- 
cal results are compared with measurements of the 3D 
motion of particles in an acoustofiuidic microchip per- 
formed by use of astigmatism particle tracking velocime- 
try (APTV) [13-15]. APTV is a very precise single- 
camera tracking method which allows a time-resolved, 
volumetric reconstruction of the trajectories of micropar- 
ticles in acoustophoretic motion. The technique is appli- 
cable to general 3D acoustophoretic motion of micropar- 
ticles influenced by both the acoustic radiation force and 
the Stokes drag from acoustic streaming. 

The paper is organized as follows. In Section II 
we derive an analytical expression of acoustic stream- 
ing in long, straight channels with rectangular cross- 
section, and we analyze the implications of this stream- 
ing for acoustophoretic motion of suspended microparti- 
cles. This is followed in Section III by a description of 
the experimental methods, in particular the astigmatism 
particle tracking velocimetry technique. In Section IV 
we compare the theoretical and experimental results for 
the acoustophoretic microparticle motion, and finally in 
Section VI we state our conclusions. 



The energy (heat), mass (continuity), and momentum 
(Navicr-Stokcs) equations take the form 

P T [d t + (v ■ V)] s = a' : Vv + V • (fc th VT) , (3a) 
d tP = -V • (pv), (3b) 

p[d t + v -V]v = -Vp+ V- »/{Vd + (Vd) t } 

+ - 1) V(*/V-w), (3c) 

where rj is the dynamic viscosity, j3 is the viscosity ra- 
tio, which has the value 1/3 for simple liquids [16], k th 
is the thermal conductivity, and a' is the viscous stress 
tensor. As in Ref. [3], we model the external ultrasound 
actuation through boundary conditions of amplitude v hc 
on the first-order velocity v 1 while keeping T constant, 



T = T 0l on all walls, 



(4a) 

v = 0, on all walls, (4b) 
n ■ v 1 = v hc (y, z) e -1 "*, added to actuated walls. (4c) 

Here n is the outward pointing surface normal vector, 
and oj is the angular frequency characterizing the har- 
monic time dependence written using complex notation. 

A. First-order fields in the bulk 



II. THEORY 

The governing perturbation equations for the thermoa- 
coustic fields are standard textbook material [16-18]. 
The full acoustic problem in a fluid, which before the 
presence of any acoustic wave is quiescent with constant 
temperature T , density p , and pressure p , is described 
by the four scalar fields temperature T, density p, pres- 
sure p, and entropy s per mass unit as well as the velocity 
vector field v. Changes in p and s are given by the two 
thermodynamic relations 



dp = jk s p dp — a p p dT, 

ds =-LdT-^dp, 
T p 



(la) 
(lb) 



which besides the specific heat capacity c p at constant 
pressure also contain the specific heat capacity ratio 7, 
the isentropic compressibility k s , and the isobaric ther- 
mal expansion coefficient a p given by 
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dp 
dp 
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Up p\dT 



(2a) 
(2b) 
(2c) 



To first order in the amplitude v hc of the imposed ul- 
trasound field we can substitute the first-order fields p Y 
and s-y in the governing equations Eq. (3) using Eq. (1). 
The heat transfer equation for T 1 , the kinematic continu- 
ity equation expressed in terms of p 1; and the dynamic 
Navier-Stokcs equation for v 1 , then become 



dtT, 




dtPi 
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Here, D th = k 



d t Pi, (5a) 
(5b) 
(5c) 

th /(/3 c p ) is the thermal diffusivity, and 



v = 7] /p Q is the kinematic viscosity. A further simpli- 
fication can be obtained when assuming that all first- 
order fields have a harmonic time dependence e~ 1LJt in- 
herited from the imposed ultrasound field Eq. (4c) . Then 
p x can be eliminated by inserting Eq. (5b), substituting 
d t p 1 = —iujpx, into Eq. (5a) and (c). Solutions of Eq. (5) 
describe the formation of thin thermoviscous boundary 
layers at rigid walls. In the viscous boundary layer of 
thickness 



(0) 



the velocity gradients are large, because the velocity field 
changes from its bulk value to zero at the walls across this 
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B. Second-order governing equations for (y 2 ) 




FIG. 1. (Color online) A cross-sectional sketch in the j/z-plane 
of the classical Rayleigh-Schlichting streaming pattern in the 
liquid-filled gap of height h between two infinite, parallel rigid 
walls (black) in the xy-plane. The bulk liquid (light shade) 
supports a horizontal standing sinusoidal pressure half-wave 
Pi (dashed lines) of wavelength A in the horizontal direction 
parallel to the walls. In the viscous boundary layers (dark 
shade) of sub-micrometer thickness 5, large shear stresses ap- 
pear, which generate the boundary-layer (Schlichting) stream- 
ing rolls (light thin lines). These result in an effective bound- 
ary condition («2y d ) (thick light arrows) with periodicity A/2 
driving the bulk (Rayleigh) streaming rolls (black thin lines). 
Only the top and bottom walls are subject to this effective 
slip boundary condition. 



In a typical experiment on microparticle acoustophore- 
sis, the microsecond timescale of the ultrasound oscilla- 
tions is not resolved. It therefore suffices to treat only the 
time-averaged equations. The time average over a full os- 
cillation period, denoted by the angled brackets (••■), of 
the second-order continuity equation and Navier-Stokcs 
equation becomes 

p Q V-(v 2 ) = -V-( Pl v 1 ), (9a) 

%V 2 (u 2 ) + /3%V(V • (v 2 )) - V(pa> 

= (Pi9*«i> +Po(( v i ■ V K) 
-<%VV)-<^iV(V- Vl )) 

- <V % • [Vv, + (V^) T ] ) 
-((f3-l)(V- Vl )V Vl ). (9b) 

Here r] 1 is the perturbation of the dynamic viscosity due 
to temperature, V = Va + Vi = v(T ) + [<9 T n(T )] T x . 
From Eq. (9) we notice that second-order temperature 
effects enter only through products of first-order fields. 
Dimensional analysis lead to a natural velocity scale u 
for second-order phenomena given by 



layer [16-18]. In water at cj/(2tt) = 2 MHz it becomes 
S « 0.4 pm. 

We focus on the transverse standing-wave resonance 
sketched in Fig. 1, which is established by tuning of u> 
in the time-harmonic boundary condition Eq. (4c) to 
achieve one of the resonance conditions n\ n /2 = w, 
n = 1, 2, 3, . . ., where A n = 27rc /w„ is the acoustic wave- 
length of the nth horizontal resonance. The associated 
first-order fields v 1 ,p 1 , and T x in the bulk of the channel 
take the form 

v i = v a sin(k n y + nn/2) e~ 1Unt e y , (7a) 
Pl = p & cos(fc„y + mr/2) e"^" 4 , (7b) 
T : =^008(^ + ^/2)6-^*, (7c) 

where k n = 1~k j\ n = mr/w is the wavenumber of the 
nth horizontal resonance, and the oscillation amplitudes 
of the first-order fields, indicated by subscript "a" , are 
related through \vjc \ - |p a /p | ~ |T a /T | <C 1, with c 
being the iscntropic speed of sound in water. The spa- 
tial form of the standing-wave resonance is determined 
entirely by the resonance frequency and the geometry of 



the resonator, while its amplitude (here v a sa 10 



[3]) 



is governed by the specific form of v hc and of the Q-factor 
of the resonance cavity. The acoustic energy density E ac 
is constant throughout the cavity and given by 



1 2 1 2 



(8) 



Po c o 



(10) 



C. The boundary condition for bulk streaming flow 

The second-order problem Eq. (9) was solved analyt- 
ically by Lord Rayleigh [19, 20] in the isothermal case 
(T = T ) for the infinite parallel-plate channel in the 
yz-plane with the imposed first-order bulk velocity v 1 , 
Eq. (7a). The resulting y-component (v 2 y ) of (y 2 ) just 
outside the boundary layers at the top and bottom walls 
becomes 



« nd > 



2y 



(11) 



as sketched in Fig. 1 for the half- wave k 1 = tt/w. In 
Rayleigh's isothermal derivation the amplitude v str of the 
streaming velocity boundary condition (v 2 y ) becomes 



3 < 
8 c n 



3 



(12) 



where the superscript "0" refers to isothermal conditions. 
Recently, Rednikov and Sadhal [21] extended this anal- 
ysis by including the oscillating thermal field T x as well 
as the temperature dependence ?7i(r) of the viscosity. 
They found that the amplitude of the streaming velocity 
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boundary condition v str then becomes 



K 1 



-K T v° stI = K T u , 
3 



7-1 



(H 

Vo a p 



th 



(13a) 
(13b) 



where the superscript " T" refers to inclusion of thermo- 
viscous effects leading to a temperature-dependent pre- 
factor multiplying the temperature-independent result. 
For water at 25 °C we find vJ tT = 1.26 v® tl . using the ma- 
terial parameter values of Table I, and in all calculations 
below we use this thermoviscous value for v str . 



D. Second-order governing equations for bulk (i>2/ 

In the bulk of the fluid the oscillating velocity and 
density fields v 1 and p x are out of phase by ir/2. Conse- 
quently (piV^ = 0, and the source term in the second- 
order continuity equation Eq. (9a) vanishes. As a result 
the time-averaged second-order velocity field (v 2 ) is di- 
vergence free or incompressible in the bulk. Hence, the 
continuity equation and the Navier-Stokes equation for 
the bulk streaming velocity field reduce to 



V-(v 2 ) = 0, 



(14a) 



(pA v i) + Po(( v i ■ V K) 
-(^V^-^iVtV-^)) 
- (V Vl ■ [Vv, + (V^) T ] ) 

-((/J-lXV-tfOVft). (14b) 



Only the y-component of the source terms on the right- 
hand side of Eq. (14b) is non-zero in the bulk, and it 
depends only on y and not on z. Consequently, their 
rotation is zero, and they can be reformulated as a gradi- 
ent term absorbed together with V(p 2 ) into an effective 
pressure gradient Vx given by, 

Vx = V(p 2 ) + (pA^i) + Po((vi ■ vh) 
-(liVV-^V-^)) 
- ■ [Vv, + (V^) T ] ) 
-((/3-l)(V-^)Vr h ). (15) 

Using this, the system of bulk equations reduces to the 
standard equation of incompressible creeping flow, 



V • (v 2 ) = 0, 



(16a) 
(16b) 



These equations together with appropriate boundary 
conditions, to be discussed below, govern the steady bulk 
streaming velocity field (v 2 ) in the microchanncl. 



TABLE I. Model parameters for water and polystyrene given 
at temperature T = 25 °C and taken from the literature as 
indicated or derived based on these. 



Parameter 


Symbol 


Value 


Unit 


Water 








Density a 


Po 


998 


kg m~ 3 


speed ot sound 


Co 


1 A OK 

1490 


m s _1 


V iscosity 


V 


U.tS93 


mPa s 


Specific heat capacity 3 


(..p 


A 1 OO 

4183 


J kg K 


Heat capacity ratio 


7 


i m a 
1.U14 




Thermal conductivity 3 


u 

" ill 


U.OUO 


W m Is. 


Thermal diffusivity 


Ah 


1.44 x 10" 7 


m 2 s- 1 


Compressibility 




448 


TPa" 1 


Thermal expansion coeff. 


Oip 


2.97 x 10" 4 


K- 1 


Thermal viscosity coeff. a 


(9tv) p 
In 


-0.024 


K- 1 


Polystyrene 








Density 13 


Pps 


1050 


kgm~ 3 


Speed of sound 


c 

"P 3 


2350 


m s _1 


Poisson's ratio d 




0.35 




Compressibility 




249 


TPa" 1 



a COMSOL Multiphysics Material Library [22], 
b Rcf. [23]. 
° Ref. [24]. 
d Ref. [25]. 



Calculated as k„ 



E. Streaming in a parallel-plate channel 

Based on Rayleigh's analysis, we first study the ana- 
lytical solution for (v 2 ) in the special case of a standing 
half wave (n = 1) in the parallel-plate channel shown 
in Fig. 1. We choose the symmetric coordinate system 
such that — w/2 < y < w/2 and — h/2 < z < h/2, and 
introduce non-dimensionalized coordinates y and z by 

2y 

y = — , with - 1< y < 1, (17a) 

w 

2z 

z = — , with - 1 < z < 1, (17b) 
h 

a = — , the aspect ratio. (17c) 
w 

In this case, using Eq. (11), the boundary conditions for 
(v 2 (y,z)) are 

M = o, 
K> = o. 

d y( v 2z) = 

where Eqs. (18c) and (18d) express the symmetry condi- 
tion at the wall-less vertical planes at y = ±1. Rayleigh 
focused his analysis of the parallel plate geometry on 
shallow channels for which a C 1. Here a — 0.4, derived 
from the aspect ratio of the microchanncl described in 
Section III and in Rcfs. [6, 8, 27], and consequently we 



sin(7ry), for 


z 


= ±i, 


(18a) 


for 


z 


= ±T 


(18b) 


for 


y 


= ±i, 


(18c) 


for 


V 


= ±i, 


(18d) 
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need to solve the case of arbitrary a. We find 

( v 2 V (y> 2 )) = v str sin(Try) (a, 5), (f 9a) 

(v2z(v,z)) =v stl .cos{iry)A- L (a,S), (19b) 

where the a- and ^-dependent amplitude functions A" 
and A 1 - for the velocity component parallel and perpen- 
dicular to the first-order wave, respectively, are given by 

A" (a, z) = S(a)| [l — 7ra coth(7ra)] cosh(7raz) 

+ naz sinh(7raz ) | , (20a) 
A 1 " (a, z) = TraB(a)^ coth(7ra) sinh(7raz) 

— z cosh(7raz) j, (20b) 

sinh(7ra) 



(a) 1 



B(a) = 



sinh(7ra) cosh(7ra) — na ' 



(20c) 



with AH(a,±l) = 1 and A x (a,±l) = 0. In Raylcigh's 
well-cited shallow-channel limit a <C 1 the amplitude 
functions reduce to 



V ' ; 2 2 



A { a ,z) k, —{z - z 



for a«l, 
* for a <C 1. 



(21a) 
(21b) 



The analytical solution of (v 2 ) for A/2 = w is illustrated 
in Fig. 2(a) and (b) for channel aspect ratios a = 1.2 
and 0.2. We note that the maximum streaming velocity 
is near the top and bottom walls. For the shallow chan- 
nel Fig. 2(b) there is furthermore a significant streaming 
velocity along the horizontal center line z = 0. However, 
the amplitude of this velocity decreases for increasing as- 
pect ratio a as shown in Fig. 2(d). 

This special case of the pure sinusoidal horizontal 
boundary condition Eq. (18a) can readily be generalized 
to any horizontal boundary condition by a Fourier ex- 
pansion in wavenumber k m = 2n/\ m = mir/w, where m 
is a positive integer, 

(V2y) = «str fit), * = ±1= (22a) 

CO 

f(y) = a m sin(m7rjf). (22b) 

TCI— 1 

As the governing equations for the second-order bulk flow 
Eq. (16) are linear, we can make a straightforward gen- 
eralization of Eq. (19), and the two velocity components 
of the superposed solution for (v 2 ^ become 

oo 

(«2j/(y> 2 )) = w str a m sm(mny) A^(ma,z), (23a) 

m— 1 
oo 

( v 2z(y,z)) = v str Y a rn C0S ( m7T y) A L (ma,z), (23b) 



>N 



(b) 

IN 

-1 

(c) 1 
in 

-1 




. - ' ' 1 - I , - ' ' ' •- . r 




FIG. 2. (Color online) Analytical results for the streaming 
velocity \v 2 ) in parallel plate channels, (a) Plot of the an- 
alytical expressions (19) and (20) for (i^) (arrows) and its 
magnitude (color plot from (black) to v BtI (white)) in the 
vertical yz cross section of a parallel-plate channel (Fig. 1) 
with A/2 = w (n = 1) and aspect ratio a = 1.2. (b) The 
same as (a), but for a = 0.2. (c) The same as (b) but for a 
standing full wave, A = w (n — 2). (d) Line plot of the am- 
plitude (v2 y (y,0)) of the streaming velocity, in units of v stI , 
along the first half of the center axis (white dashed lines in 
(a) and (b)) with A/2 = w for aspect ratios a = 0.2, 0.5, 0.8, 
and 1.2. (e) Line plot of the maximum (v2 y (y,0)) rnax of the 
center-axis streaming velocity, in units of w str , as function of 
aspect ratio for the resonances nA/2 = w, with n = 1, 2, and 
3, respectively. 
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where the wave index m multiplies both the horizon- 
tal coordinate y and the aspect ratio a. Note that 
A" (ma, ±1) = 1 and A x (ma,±l) = 0. The result- 
ing steady effective pressure x is j us t the weighted sum 
of the partial pressures x m °f each Fourier component, 

m— 1 ^mXm ' 

In Fig. 2(c) is shown the streaming velocity field for 
the higher harmonic boundary condition f(y) = sin(mn/) 
with n = 2. Furthermore, Fig. 2(e) shows how the maxi- 
mum (v 2y (y, 0)) max of the center-axis streaming velocity 
decays as function of aspect ratio a for n = 1, 2, and 3. 
Given sufficient room, the flow rolls decay in the verti- 
cal direction on the length scale of A„/4. Since n is the 
number of half wavelengths of the first-order resonance 
pressure across the channel, we conclude that the stream- 
ing amplitude in the center of the channel decreases for 
higher harmonic resonances. 



F. Streaming in a rectangular channel 

Moving on to the rectangular channel cross section, we 
note that the only change in the problem formulation is to 
substitute the symmetry boundary conditions Eqs. (18c) 
and (18d) by no-slip boundary conditions, while keep- 
ing the top-bottom slip boundary conditions Eqs. (18a) 
and (18b) unaltered, 







sin(Try), for z 


= ±1, 


(24a) 


( V 2z) 


= o, 


for z 


= ±1, 


(24b) 




= 0, 


for y 


= ±1, 


(24c) 




= 0, 


for y 


= ±1. 


(24d) 



If we want to use the solution obtained for the parallel- 
plate channel, we need to cancel the vertical velocity com- 
ponent (v 2z ) on the vertical walls at y = ±1. This leads 
us to consider the problem rotated 90°, where the first- 
order velocity field is parallel to the vertical walls (in- 
terchanging y and z), and the fundamental wavelength 
is A/2 = h, and the aspect ratio is w/h = a -1 . As the 
governing equations for the bulk flow Eq. (16) are linear, 
we simply add this kind of solution to the former so- 
lution and determine the Fourier expansion coefficients 
such that the boundary conditions Eq. (24) are fulfilled. 
Given this, Eq. (23) generalizes to 



a m sin(m7ry) (ma, z) 



+ b m A (ma ,y) cos(rrmz) 



a m cos(m7ry) A ± (ma, z) 



(25a) 



m— 1 

+ 6 m A"(ma _1 ,y) sin(ra7T2) 



(25b) 



The two pcrpcndicular-to-the-wall velocity conditions 
Eqs. (24b) and (24c) are automatically fulfilled as they 
by construction are inherited from the original conditions 
Eqs. (18b) and (18c). The unknown coefficients a m and 
b m are thus to be determined by the parallel-to-the-wall 
conditions Eqs. (24a) and (24d). 

Using (v 2y ) in the form of Eq. (25a), boundary condi- 
tion Eq. (24a) becomes 



sin(Try) = ^ 



a m sm(miry) + (-l) m b m A J -(ma 1 ,y) 



(26) 

The discrete Fourier transform of this equation, i.e. mul- 
tiplying by sin(j7ry), where j is an arbitrary integer, and 
integrating over y from —1 to 1, becomes 

oo 

hi = E [ 5 i,m a ™+ A ^( a_1 ) 6 ™] , 3 = 1, 2, 3, . . . (27) 

m— 1 

where the (j, m)'th element Af m of the a-dependent ma- 
trix A 1 ^ is given by 



A£ m (0 = (-l) m / dyA^(ma-\y) sin(jny). 



(28) 



Introducing the coefficient vectors a and b and the first 
unit vector e 1 with m'th components a m , b m , and 8 1 m , 
respectively, we can write Eq. (27) as the matrix equation 



e, = a 



A ± (a" 1 ) b. 



(29) 



Likewise, using Eq. (25b) and multiplying it by sin(jnz), 
where j is an arbitrary integer, and integrating over z 
from —1 to 1, the zero-parallel-component boundary con- 
dition Eq. (24d) can be written as the matrix equation 



= A ± (a) ■ a + b. 



(30) 



Solving the equation system Eqs. (29) and (30), the co- 
efficient vectors a and b becomes 



a= I- A ± (a- 1 )A ± (a) 
b = -A ± (a) ■ a. 



-i ■ 



(31a) 
(31b) 



A comparison between results for the classical parallel- 
plate geometry and the new results for the rectangular 
geometry is shown in Fig. 3. It is seen how the velocity 
profile of the rectangular channel solution, Eq. (25), is 
suppressed close to the wall in comparison to the parallel- 
plate channel solution, Eq. (19). Note that for the nth 
resonance, k n = mr/w, the unit vector e 1 in Eq. (31a) is 
replaced by (— l) n_1 e n , with the sign originating from 
the n-dependent phase shift in the streaming boundary 
condition Eq. (11). 
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FIG. 3. (Color online) Analytical results comparing the 
streaming velocity field (v 2 ) in the parallel plate and the rect- 
angular channel, (a) Color plot from (black) to v stI (white) 
of the analytical expression for (« 2 ) Eqs. (19) and (20) in 
the classical parallel-plate geometry with a half-wave reso- 
nance A/2 — w (n = 1). Due to symmetry, only the left half 
( — 1 < y < 0) of the vertical channel cross section is shown. 

(b) As in (a) but for (u 2 ) in the rectangular channel Eqs. (25) 
and (31), including the first 20 terms of the Fourier series, (c) 
Line plots of {v2 y (y,0)') in units of v stI along the left half of 
the center line for the parallel-plate channel (dashed lines) and 
the rectangular channel (full lines) for aspect ratios a — 0.1, 
0.4, and 0.8 and the half-wave resonance A/2 = w. (d) As in 

(c) but for the full- wave resonance A = w (n — 2). 



G. Acoustophoretic particle velocity 

The forces of acoustic origin acting on a single mi- 
croparticlc of radius a, density p p , and compressibility 
K p undergoing acoustophoresis with velocity it p in a liq- 
uid of density p , compressibility k s , and viscosity rj , are 
the Stokes drag force _F dra s = 67r?/ a[(i> 2 ) — it p ] from the 
acoustic streaming (v 2 ) and the acoustic radiation force 
F* . Given an observed maximum acoustophoretic ve- 
locity of u p < 1 mm/s for the largest particles of di- 
ameter 2a = 5.0 pm, the Reynolds number for the flow 
around the particle becomes p Q 2au p /rj < 6 x 10~ 3 , and 
the time scale for acceleration of the particle becomes 
Tacc = [(4/3)7ra 3 /9 p ] / [677770] w 2 us. Since the acceler- 
ation time is much smaller than the time scale for the 
translation of the particles r trans = w/(2u p ) > 0.1 s, the 
inertia of the particle can be neglected, and the quasi 
steady-state equation of motion, _F dra s = —F Tad } for a 



spherical particle of velocity tt p then becomes 

^rad 



67T?7 a 



+ (v 2 ) = u™ d + (v 2 ) 



(32) 



where it rad is the contribution to the particle velocity 
from the acoustic radiation force. The streaming velocity 
(w 2 ) is given in the previous subsections, while an analyt- 
ical expression for the viscosity-dependent timc-avcraged 
radiation force F rad in the experimentally relevant limit 
of the wavelength A being much larger than both the 
particle radius a and the boundary layer thickness 5 was 
given recently by Settnes and Bruus [28] . For the case of a 
ID transverse pressure resonance, Eq. (7b), the viscosity- 
dependent acoustic radiation force on a particle reduces 
to the x- and z-independent expression 

F rad (y) = ina%E ac $(k,p,~6) sin [nir(y + l)]e tf . (33) 

The acoustic contrast factor $ is given in terms of the 
material parameters as 



= 1 - «, 

2[i-r(5)](p-i) 



/ 2 GM) = 

T(S) = 



2p 
3 

~2 



+ 1 - 3T(S) 
l + i(l + <5) 



(34a) 
(34b) 

(34c) 
(34d) 



where k = K p /n s , p = p p /p , and 5 = 5/a. Using 
Eq. (33) for the transverse resonance, u rad only has a 
horizontal component u 



rad 



rad 



u sin [mr(y + 1)] , n = 1, 2, 3, . 



(35a) 



where the characteristic particle radius a is given by 




(35b) 



with S given by Eq. (6). The acoustophoretic particle ve- 
locity u p will in general have a non-zero z-componcnts, 
due to the contribution from the acoustic streaming (f 2 ) ■ 
However, for the special case of particles in the horizon- 
tal center-plane z = of a parallel-plate or rectangular 
channel, the vertical streaming velocity component van- 
ishes, (v2z(y,0)) = 0. From Eqs. (19a) and (35a) we find 
that the horizontal particle velocity component u p (y,0) 
in a parallel-plate channel is given by the sinusoidal ex- 
pression, 



uUy,0) = u Ci 



^j-K T AHna,0) 



sm 



[»7r(jf + l)]. (36) 



Since by Eq. (20a) A" (na, 0) is always negative, it follows 
that the streaming-induced drag and the radiation force 
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have the same direction in the horizontal center plane of 
the channel. For the rectangular channel using Eq. (25a), 
the expression for u%(y,0) becomes 



u%(y,o) = u o\ -j sin [ nn (y + !)] 



(37) 



m— 1 



a m sin(m7ry) (ma, 0) + b m A ± (ma 1 , y) 



which is not sinusoidal in y but still proportional to u . 
This particular motion in the ultrasound symmetry plane 
is studied in detail in Ref. [27]. 



III. EXPERIMENTS 

We have validated experimentally the analytical 
expressions derived above by measuring trajectories 
of micrometer-sized polystyrene particles displaced by 
acoustophoresis in a long, straight silicon/glass mi- 
crochannel with rectangular cross section. A fully three- 
dimensional evaluation of the particle trajectories and 
velocities was performed by means of the astigmatism 
particle tracking velocimetry (APTV) technique [13, 14] 
coupled to the temperature-controlled and automated 
setup presented in Rcf. [8]. APTV is a single-camera 
particle-tracking method in which an astigmatic aber- 
ration is introduced in the optical system by means of a 
cylindrical lens placed in front of the camera sensor. Con- 
sequently, a spherical particle image obtained in such a 
system shows a characteristic elliptical shape unequivo- 
cally related to its depth-position z. More details about 
calibration and uncertainty of this technique, as well as 
comparison with other whole-field velocimetry methods 
for microflows, can be found in Refs. [14, 15]. 



A. Acoustophoresis microchip 

The acoustophoresis microchip used for the experi- 
ment was the one previously presented in Refs. [6, 8, 27]. 
Briefly, a rectangular cross section channel (L — 35 mm, 
to = 377 pm, and h = 157 pm) was etched in silicon. A 
Pyrex lid was anodically bonded to seal the channel and 
provided the optical access for the microscope. The outer 
dimensions of the chip are L = 35 mm, W = 2.52 mm, 
and H = 1.48 mm. From top and down, glued together, 
the chip was placed on top of a piezoceramic transducer 
(piezo) , an aluminum slab to distribute heat evenly along 
the piezo, and a Peltier element to enable temperature 
control based on readings from a temperature sensor 
placed near the chip on the transducer. This chip-stack 
was mounted on a computer-controlled xyz-st&ge. Ultra- 
sound vibrations propagating in the microchip were gen- 
erated in the piezo by applying an amplified sinusoidal 
voltage from a function generator, and the resulting piezo 
voltage 17 was monitored using an oscilloscope. 



B. APTV set-up and method 



The images of the particles in the microfluidic chip 
were taken using an epifluorescent microscope (DM2500 
M, Leica Microsystems CMS GmbH, Wetzlar, Germany) 
in combination with a 12-bit, 1376x1040 pixels, inter- 
line transfer CCD camera (Sensicam QE, PCO GmbH). 
The optical arrangement consisted of a principal ob- 
jective lens with 20 x magnification and 0.4 numerical 
aperture and a cylindrical lens with focal length f cy \ = 
150 mm placed in front of the CCD sensor of the cam- 
era. This configuration provided a measurement volume 
of 900 x 600 x 120 pm 3 with an estimated uncertainty 
in the particle position determination of ±1 pm in the in- 
direction and less than ±0.1 pm in the x- and y-dircction. 
Two scan positions along the z-direction were used to 
cover the whole cross-sectional area of the channel. 

Monodisperse spherical polystyrene particles with 
nominal diameters of 5.33 pm (SD 0.09) and 0.537 pm 
(PDI 0.005) were used for the experiments (p ps = 
1050 kg m~ 3 and K ps = 249 TPa" 1 ). For simplicity we 
will refer to them as 5-pm-diam and 0.5-pm-diam par- 
ticles, respectively. The particles were fabricated and 
labeled with a proprietary fluorescent dye by Micropar- 
ticles GmbH to be visualized with an epifluorescent mi- 
croscopy system. The illumination was provided by a 
continuous diode-pumped laser with 2 W at 532 nm wave- 
length (www.mylaserpage.de). 

Once the particle 3D positions had been detected us- 
ing the APTV technique, their trajectories and velocities 
were calculated. Due to the low seeding density in the 
experiments, the particle inter-distance was large enough 
to employ a simple nearest-neighbor approach in which 
the particle in one frame is identified with the closest par- 
ticle in the next frame. The method was compared with 
more sophisticated ones as predictors and probabilistic 
algorithms with identical results. Trajectories composed 
by less than 5 particle positions were rejected. From 
the obtained trajectories the velocities could be calcu- 
lated given the capture rate of the camera. Different 
approaches have been followed depending on the type 
of trajectories expected. For particles following almost 
straight paths as those dominated by radiation force, a 
simple 2-position approach was used and the velocities 
were calculated based on the frame-to-frame particle dis- 
placement. For particles with highly curved paths, like 
those present in streaming-dominated flows, a more so- 
phisticated multi-frame approach has been followed, as 
those reported already by Hain and Kahler [29] for pPIV. 
In our case, each velocity data point was calculated from 
a trajectory segment composed by 4 consecutive points. 
Such multi-frame approach applied for PTV has been 
shown to better solve the velocity vector positions and 
values when the trajectories present large curvatures and 
for high-shear flows [30]. 
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IV. RESULTS 
A. APTV measurements 

Examples of the measured 3D trajectories of the 5- 
pm-diam particles are shown in Fig. 4(a). The data 
was collected from 10 consecutive experiments with the 
piezo operated at 1.94 MHz and peak-to-peak voltage of 
U pp = 0.91 V. An overall number of 111 trajectories were 
determined. The 5-pm-diam particles are affected mainly 
by the acoustic radiation force that quickly pushes 
them to the center of the channel with a horizontal veloc- 
ity «p [4, 27]. At the vertical pressure nodal plane y = 0, 
F d vanishes and the hitherto negligible drag force from 
the acoustic streaming, shown in Fig. 2(b), slowly drags 
the particles towards the top and bottom of the channel. 

Examples of the measured 3D trajectories of the 0.5- 
um-diam particles are shown in Fig. 4(b). The data 
was collected from four consecutive experiments with the 
piezo operated at 1.94 MHz and pcak-to-pcak voltage of 
U = 1.62 V. An overall number of 731 trajectories were 
determined. The acoustic radiation force F^ ad is in this 
case minute and the particles are primarily transported 
by the acoustic streaming (v 2 ) of the fluid resulting in 
particle trajectories following the four vertical vortices in 
the bulk, shown in Fig. 2(b). 



B. Comparison of theory and experiments 

Theoretically, the acoustophoretic particle velocity it p 
is given by Eq. (32) combined with the expressions for 
the streaming velocity of the liquid Eqs. (25) and (31) 
and the expression for the radiation force on the particles 
Eq. (33). The amplitudes of both the acoustic streaming 
and the radiation force depend linearly on the acoustic 
energy density E ac through Eqs. (12) and (33). To make 
a theoretical prediction of the motion of the 0. 5-pm-diam 
particles we need to determine the acoustic energy den- 
sity E® c 5 11111 . This calibration is done in situ based on 
the measurements of the 5-pm-diam particles, by the fol- 
lowing three-step procedure. 

First, we determine the acoustic energy density E ac mi 
for the experiment with the 5-um-diam particles. This is 
done by fitting the sin(7ry)-dependent expression Eq. (36) 
for u^(y, 0) to the measured instantaneous velocities, us- 
ing the amplitude as the only fitting parameter [6, 27]. 
The small contribution from the acoustic streaming to 
the 5-um-diam-particle velocity is taken into account al- 
though it constitutes only 6 % of the total particle veloc- 
ity. The fit showed good agreement between theory and 
experiment, and after correcting for a wall-enhanced drag 
coefficient of 1.032 at the horizontal center plane (see 
Refs. [7, 27, 31, 32]), we extracted the acoustic energy 
density E^ lm = (20.6 ± 0.7) J/m 3 , where the la stan- 
dard error of the estimated value is stated. Since both 
the wall-enhanced drag coefficient and the drag force 




V [Mm] 200 

FIG. 4. (Color online) Measured particle trajectories (thin 
black lines) obtained using the 3D-APTV technique in the 
microchannel (gray walls) actuated at the 1.94-MHz horizon- 
tal half-wave resonance. For selected trajectories, the parti- 
cle positions are represented by dots, (a) 5-um-diam parti- 
cles moving (red arrows) to the vertical center plane y = 0, 
and (b) 0.5-um-diam particles exhibiting circular motion as 
in Fig. 2(b). 



from the acoustic streaming fluid velocity are height- 
dependent, only five trajectories of 5-pm-diam particles 
close to the horizontal center line (z = 0) qualified for 
use in the fit, based on a criterion of |z | < 6 pm. 
The starting positions (x ,y , z ) of the five tracks were 
(34 pm, — 115 pm,6 pm), (310 pm, — 66 pm, — 6 pm), 
(482 pm, — 35 pm, — 5 pm), (74 pm, 115 pm, 2 pm), and 
(350 pm, 128 pm, pm), and they all reached the vertical 
center plane y = 0. 

Second, the acoustic energy density E®f 11111 for the 
experiment with the 0. 5-pm-diam particles was deter- 
mined, using the result for -E^c 1 ™ combined with the fact 
that E ac scales as the square of the applied voltage U 
[6]. The measured voltages for the two experiments are 
U°p 11111 = (1.62 ± 0.01) V and [/pp 1 ™ = (0.91 ± 0.01) V, 
where the stated error corresponds to the standard de- 
viation of a series of voltage measurements, with the 
power turned off in between each measurement. The 
derived value for E^ 11111 , taking into account the er- 
rors of J7pp 5 11111 , Up P im \ and E^ m , becomes E° c 5 11111 = 
(Upp 5 11111 /C/pp 11111 ) 2 El 11111 = (65 ± 2) J/m 3 , with the con- 
tribution from the error of the measured voltages being 
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FIG. 5. (Color online) Comparison between experimental, 
analytical, and numerical studies of the acoustophoretic par- 
ticle velocities u p of 0.5-um-diam polystyrene particles in wa- 
ter. The particle velocities it p (vectors) and their magnitude 
(color plot ranging from um/s (black) to 63 um/s (white) in 
all three plots), are shown in the vertical cross-section of the 
microchannel, divided into a pixel array consisting of 37-by-15 
square bins of side length 10 p. The axes of the plot coincide 
with the position of the channel walls, (a) The APTV mea- 
surements of the 0.5-um-diam particles, shown in Fig. 4(b), 
projected onto the vertical cross section. The maximum ve- 
locity is 63 um/s. Close to the side walls experimental data 
could not be obtained, which is represented by hatched bins, 
(b) Analytical prediction of u p based on Eq. (32), taking both 
the radiation force and the streaming-induced drag force into 
account. The first 20 terms of the Fourier series for (v 2 ), 
Eq. (25), have been included in the calculation. The maxi- 
mum velocity is 59 um/s. There are no free parameters in 
this prediction as the acoustic energy density was calibrated 
in situ based on measurements of large 5-um-diam particles, 
shown Fig. 4(a). (c) Numerical validation of the analytical 
result for tt p using the method described in Muller et al. [3]. 
The numerical solution has been scaled by the thermoviscous 
pre- factor to the streaming amplitude Eq. (13). The maxi- 
mum velocity is 59 um/s. 



negligible. 

Third, based on Eq. (10), the derived value for the en- 
ergy density E® c 5 11111 is used in the analytical expression 
for the particle velocities Eq. (32). The radiation force 
is given by Eq. (33) and the acoustic streaming velocity 
is given by Eqs. (25) and (31), using the thcrmoviscous- 




FIG. 6. (Color online) (a) Color plot of the number of times 
the velocity has been measured in each square bin. (b) Color 
plot of standard error of the mean (SEM) particle velocity in 
each square bin. 



corrected amplitude Eq. (13). The contribution from 
the acoustic radiation force to the 0.5-um-diam-particle 
velocity is small and constitutes only 12 % of the total 
particle velocity in the horizontal center plane z = 0. The 
contribution from the radiation force to the 0.5-um-diam- 
particlc velocity is not corrected for the wall-enhanced 
drag coefficient, since this is minute for these small par- 
ticles. 

To compare the experimental results and the analytical 
prediction, we consider the 0.5-um-diam-particle veloci- 
ties in the vertical cross section, yz-planc, of the channel 
as in Fig. 1, 2, and 3. In Fig. 5 are shown color plots 
of (a) the experimentally measured acoustophoretic ve- 
locities for the 0.5-um-diam particles, (b) the analytical 
prediction of the same, and (c) the numerical validation 
of the analytical result using the methods of Muller et al. 
[3] . The three data sets are shown on the same 37 x 15 bin 
array and with the same color scale. The experimental 
and the analytical velocities agree well both qualitatively 
and quantitatively, although the experimental velocities 
are approximately 20 % higher on average. The exper- 
imental results for the particle velocities, Fig. 5(a), is 
found as the mean of several measurements of the par- 
ticle velocity in each bin. The number of measurements 
performed in each bin is shown in Fig. 6(a), while the 
standard error of the mean (SEM) particle velocity is 
shown in Fig. 6(b). These plots show that we typically 
have between 20 and 70 measurements in each bin and 
the experimental error is on average only 1 um/s, while 
the relative experimental error is on average 4 %. The 
error of the theoretical prediction is given by the relative 
error of 4 % on the estimated value for the energy density 
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FIG. 7. (Color online) (a) Color plot of the difference be- 
tween the experimental and analytical acoustophoretic par- 
ticle speeds, Au p Eq. (38). (b) Line plots of A« p along the 
dashed lines in (a), marked A, B, C, D, E, and F, with error 
bars indicating the la error of Au p . The lines are positioned 
at y — pm, y = ±91.7 p, 2 = p, and z = ±52.3 um. 
The off-center lines go through the rotation centers of the 
flow rolls, and consequently it p ~ u % e y m B, D, and F, while 
u p w u P e z in A, C, and E. 

The quantitative differences between the experimen- 
tal particle velocities Fig. 5(a) and the analytical pre- 
diction Fig. 5(b) are emphasized in Fig. 7, showing the 
difference Aw p between the experimental and analytical 
acoustophoretic particle speeds 

A^ = K p |-|< nl |. (38) 

We have chosen to consider the difference of the abso- 
lute velocity values, |« p xp | — l^anil' instead of the abso- 
lute of the difference, \u p xp — u^ nl \, because the former 
allows us to see when the experimental velocity respec- 
tively overshoots and undershoots the analytical predic- 
tion. Fig. 7(a) shows a color plot of Au p in the channel 
cross section, while Fig. 7(b) shows line plots of Au p 
along the dashed lines in Fig. 7(a), allowing for more de- 
tailed study of the spatial dependence of the difference. 
These lines arc chosen to go through the rotation centers 
of the flow rolls. The error bars in Fig. 7(b) show the 
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FIG. 8. (Color online) Experimental data from Ref. [27] com- 
pared with the theoretical predictions of Eqs. (36) and (37). 
uPIV measurements, in the center plane z = 0, of the y- 
component of the acoustophoretic velocity (u p (y,Q)} x (open 
and closed dots) for 0.6-um-diam polystyrene particles in wa- 
ter, small enough that streaming dominates and it p m (v 2 )- 
The observed motion (red arrows) in (a) and (b) resembles the 
analytical results shown in Fig. 2(b) and (c), respectively. For 
each value of y, the measured velocity w p is averaged along 
the K-coordinate, with resulting SEM smaller than the size of 
the dots. The sinusoidal parallel- plate prediction (thin line), 
Eq. (36), is fitted to the data points far from the side walls 
(open dots), while the rectangular-channel prediction (thick 
line), Eq. (37), is fitted to all data points (open and closed 
dots). In both fit the acoustic energy, _E ac , is treated as a free 
parameter, (a) The half- wave resonance A/2 = to (n = 1) 
with / = 1.940 MHz and U pp = 1 V. (b) The full-wave reso- 
nance A = w (n = 2) with / = 3.900 MHz and U pp = 1 V. 



lcr error of Ait p , taking into account both the SEM for 
the experimental measurements, Fig. 6(b), and the error 
of the analytical prediction (4 %) inherited from the de- 
rived value for E® c 5 11111 . The experimental and analytical 
velocities do not agree within the error of Au p , moreover, 
a trend of the experimental velocities being larger than 
the analytical predictions is seen. 

A further comparison between the analytical model 
presented in this paper and experimental measurements 
on 0.6-um-diam polystyrene particles from Ref. [27] is 
shown in Fig. 8. These particles are dominated by the 
drag from the acoustic streaming, and in this compari- 
son we are only interested in studying how the side walls 
influence the shape of u p (y,Q), Eq. (37). Consequently, 
the amplitude of the streaming velocity, and thus the 
acoustic energy density, is treated as a fitting parame- 
ter. The experimental results support our analytical pre- 
diction Eq. (37) (thick line) for the rectangular channel 
with side walls, which shows a suppression of u p near the 
walls compared to the sinusoidal form of u p in Eq. (36) 
(thin line) predicted for the parallel-plate channel with- 
out side walls. This is particularly clear for the full wave 
resonance A = w (n — 2) Fig. 8(b). The difference in 
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the amplitude of vP between Fig. 8(a) and (b) is due to 
differences in the resonance modes, i.e. E &c is not the 
same even though [/ is. 



V. DISCUSSION 

The comparison shows good agreement between the 
experimental measurements and the analytical prediction 
of the streaming-induced particle velocities. The qualita- 
tive agreement is seen in Fig. 5 for the two-dimensional 
topology of the particle motion, and in Fig. 8 for the 
effect of the side walls. Quantitatively, the experimen- 
tal and analytical results agree within a mean relative 
difference of approximately 20 %, a low deviation given 
state-of-the-art in the field. However, as illustrated by 
the statistical analysis in Fig. 7, the differences Au p are 
larger than the estimated lcr-errors. This could indicate 
a minor systematic error in the experimental procedure 
or in the theoretical model, or be due to underestima- 
tion of the experimental error involved in the analytical 
prediction. 

In the 5-um-diam-particle experiment the acoustic en- 
ergy density is determined using only five particle tra- 
jectories close to the channel center z = 0. This is 
reasonable as the error of the calculated energy den- 
sity is relatively low (3%), however, a calculation based 
on more particle trajectories would be desirable. This 
can be realized through more experimental repetitions 
or through implementation of the 2D-dependence of the 
wall-enhanced drag force, allowing for use of off-center 
particle trajectories. One source of error that has not 
yet been discussed arises from the assumption made that 
the acoustic energy density, and thus the acoustic forces, 
does not depend on the x-position in the investigated field 
of view. In the same setup, Augustsson et al. [8] observed 
negligible field gradients in the x-direction in some field of 
views and significant ones in others. This inhomogeneity 
was considered here when making the measurements: we 
made sure to check that the five 5-um-diam-particle tra- 
jectories sample the x-range reasonably well and exhibit 
only negligible variations in the acoustic energy density 
function of x-position. 

In the 0. 5-um-diam-particle experiment the statistics 
and sampling of the x-range are good, however they could 
still be improved to achieve better statistics close to the 
walls. The relative positions of the 0.5-um-diam particles 
are accurately determined by use of the APTV technique, 
whereas the absolute position in the channel, which was 
used to compare with theory, is difficult to determine 
precisely and might also be improved. Furthermore, ac- 
curate measurements of the channel dimensions are also 
important, as these are key parameters in the theoretical 
model. 

The analytical model could be improved in several 
ways. The treatment of the liquid could be extended 
by including thermal dependence of more material pa- 
rameters such as the specific heat capacity ratio 7, ther- 



mal expansion a p , compressibility n s , and speed of sound 
c . The influence of the surrounding chip material could 
be included, thereby relaxing the assumptions of infinite 
acoustic impedance (ideal reflection) and infinite thermal 
conduction (ideal heat sink) of the channel walls. Solving 
the full elastic wave problem in the whole chip is beyond 
analytical solutions, but is, however, possible with nu- 
merical models. This might be necessary to achieve ac- 
curate quantitative agreement between theoretical pre- 
dictions and experiments. Furthermore, the analytical 
and numerical models assume an ideal rectangular chan- 
nel cross section, which is crucial since the generating 
mechanism for the acoustic streaming takes place within 
the um-thin acoustic boundary layer. Even small defects, 
such as uneven surfaces on the um-scale, might lead to 
changes in the acoustic streaming velocity field. 



VI. CONCLUSIONS 

In this work we have for a rectangular microchannel 
derived an analytical expression for the acoustophoretic 
velocity of microparticles resulting from the acoustic ra- 
diation force and the acoustic streaming-induced drag 
force, and wc have successfully compared it with a direct 
numerical solution of the governing equations. We have 
also accurately measured 3D trajectories of 0.5-um-diam 
and 5-um-diam particles in an acoustically actuated mi- 
crochannel, with an average relative experimental error 
of 4 % for the 0.5-um-particle velocities. This allowed us 
to perform a quantitative comparison in 3D between the- 
ory and experiments of streaming-induced particle veloc- 
ities in a rectangular channel. The analytical derivation 
successfully predicted the measured streaming-induced 
0. 5-um-diam-particle velocities, with qualitative agree- 
ment and quantitative differences around 20 %, a low 
deviation given state-of-the-art in the field. This shows 
that the time-averaged second-order perturbation model 
of the governing equations yields an adequate description 
of the acoustophoretic particle motion. 

The differences between the theoretical prediction and 
the experimental results emphasize the need for further 
extensions of the analytical model, along with improved 
numerical simulations [3]. Aiming for more detailed 
quantitative studies of acoustophorcsis, the results also 
stress the need for improved accuracy of the measure- 
ments of the channel dimensions and the absolute posi- 
tions of the particles in the microchannel. The trinity 
of analytical, numerical, and experimental studies of the 
acoustophoretic particle motion enhance the understand- 
ing of acoustophoresis and supports a more elaborate and 
broader application of acoustophoresis. 
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